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Abstract 



We perform the complete bosonization of 2+1 dimensional QED with one 
fermionic flavor in the Hamiltonian formalism. The Fermi operators are explic- 
itly constructed in terms of the vector potential and the electric field. We carefully 
specify the regularization procedure involved in the definition of these operators, 
and calculate the fermionic bilinears and the energy - momentum tensor. The 
algebra of bilinears exhibits the Schwinger terms which also appear in perturba- 
tion theory. The bosonic Hamiltonian is a local, polynomial functional of Ai and 
Ei, and we check explicitly the Lorentz invariance of the resulting bosonic the- 
ory. Our construction is conceptually very similar to Mandelstam's construction 
in 1+1 dimensions, and is dissimilar from the recent bosonization attempts in 2+1 
dimensions, which hinge crucially on the presence of a Chern - Simons term. 
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1 Introduction 



It is well known that, in 1+1 dimensions, fermion field operators can be constructed in 
terms of local bosonic fields. This bosonization procedure is of great theoretical interest 
in itself, since it provides a method of mutually mapping different local quantum field 
theories with a priori different Hilbert spaces. In the pioneering works of Coleman and 
Mandelstam 0, the bosonization of a single Dirac field was performed. Subsequently it 
was realized, that this procedure is much more general. An extension was found to the- 
ories of several Fermi fields 0, and eventually non-Abelian bosonization was introduced 
by Witten 0]. Bosonization has proved to be very helpful in the analysis of different 
1+1 dimensional models. It has been used in a variety of contexts: relativistic quantum 
field theory, condensed matter systems, string theory. 

From the point of view of particle physics, bosonization techniques have found per- 
haps their most interesting applications in the realm of gauge theories. There they 
provide insights into the strong coupling regime of both Abelian |^, and non-Abelian ^ 
gauge theories, as well as better understanding of t'Hooft's solution of large N QCD2 
0. Also in many other instances, in which 1+1 dimensional theories are used to model 
the behaviour of QCD4, bosonization has proven helpful |Q. 

Clearly, the extension of bosonization to higher dimensions is highly desirable. It 
is reasonable to hope that, apart from its intrinsic theoretical interest, bosonization 
in higher dimensions will also find diverse applications. Although it is not likely that 
bosonized theories in 2+1 and 3+1 dimensions will be exactly solvable, they may be 
amenable to different analytical methods than their original fermionic formulations. A 
local bosonic formulation of gauge theories with fermions should also be very interesting 
from the lattice gauge theory point of view, since it would provide a convenient starting 
point for numerical simulations. 

In the late 1970's and early 1980's several attempts, using different approaches, have 
been made to perform bosonization in more than 1+1 dimensions. None of them however 
has lead to a local bosonic theory. One approach was to perform an exact analog of the 
Jordan - Wigner transformation for lattice fermions |Q . The resulting bosonic theory is 
a Z2 gauge theory, and the bosonic variables are not gauge invariant, and consequently 
nonlocal. The nature of the continuum limit of this bosonic model is also not known. In 
another approach [|l^, one divides a two (or three) dimensional space into one dimen- 
sional subspaces, and performs one - dimensional bosonization in each subspace. This 
procedure, however, does not give local expressions for some of the local fermionic bilin- 
ears. In particular, the fermion mass term is given by a complicated nonlocal operator. 
Related attempts, using tomografic projection |jll|, have produced similar results. 

After the relative failure of these early approaches, it has been tacitly assumed untill 
quite recently, that the extension of bosonization to higher dimensions is impossible. 
The feeling was that 1+1 dimensions is a very special case, since there is no spin in 1+1 
dimensions, and therefore no "real" difference between bosons and fermions. 

The renewed interest in this problem was triggered by a possible connection between 
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the phenomenon of Fermi - Bose transmutation in 2+1 dimensions, and novel condensed 
matter systems, most notably high Tc superconductors. Polyakov argued |12], that 
the addition of a Chern - Simons term [|I3] to the Lagrangian of QED3 changes the 
statistics of charged excitations. His argument was elaborated further in many works 
13, and the explicit realization of this Chern - Simons induced mechanism of Fermi - 



Bose transmutation in lattice theories has been given in |jT5|, [0. 

There is, however, an important conceptual difference between the Mandelstam con- 
struction in 1+1 dimensions, and the Chern - Simons bosonization in 2+1 dimensions. 
Let us recall the basic bosonization formulae in 1+1 dimensions [Q. The fermionic op- 
erator ifja is constructed in terms of a scalar bose field 0, and its conjugate momentum 



as follows, 



^'1,2(3;) = exp 



-2i 



vr 



dyp{y) ± 



which gives the following expressions for the bilinears, 

^TmV" = --^(^t^udu(j); i>iJ oc: cos/ 
Zn 



(1) 



(2) 



One important aspect of these formulae is that the fermion number current, which in the 
original theory is conserved due to the equations of motion, bosonizes into a topological 
current, which is conserved trivially. The U(l) charge, when expressed in terms of 0, is 
a topological, rather than a Noether, charge. The bosonization is therefore achieved by 
a duality transformation: one constructs the fermionic operators ip which carry a global 
U(l) fermion number charge, in terms of a local bosonic field (p, which is itself neutral. 



The 2+1 dimensional construction of |15 
however 



16| is very different in this respect. (See, 



). In that case the expression for the Fermi field is 



ip{x) = (f){x)exp 



i / d^y9{x - y)p{y) 



(3) 



where 9{x) is the planar angle, and p is the charge density operator. The fermi field ip 
and the bose field carry the same global quantum numbers. The element of duality is 
notably missing here. Moreover, in this picture, a free fermion is bosonized into a particle 
interacting with the vector potential with Chern - Simons action. But the associated 
bose field is not gauge invariant, and, after gauge fixing, becomes nonlocal. So one does 
not really construct fermionic operators in terms of local bosonic fields. 

This 2+1 dimensional construction hinges crucially on the presence of the Chern - 
Simons term in the action. The Fermi - Bose transmutation occurs only for a fixed value 
of its coefficient. For other values of this coefficient this procedure gives nonlocal anyonic 
fields with fractional statistics. Another unsatisfactory feature of the procedure is, that 
it has been consistently implemented only for nonrelativistic fermions [|T^, |T^. When 
applied to continuum relativistically invariant theories |jl8|, it suffers from regularization 



ambiguities, and it is not clear how to interpret the formal results. For example, since 
[p(?/), 0(a;)] = (j){x)5'^{x — y), and the angle function is not defined at the origin, the 
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ordering on the right hand side of expression @ is completely ambiguous, and formal 
manipulations using it are ill- defined. Furthermore, the covariant Dirac field has not 
yet been constructed in this framework, even formally. 

In this paper we take a different approach to bosonization. Our aim is to construct a 
Dirac doublet of fermionic operators in 2+1 dimensions in terms of a local bosonic field. 
The main hypothesis on which we base this construction is the following. We assume 
that, in terms of the bose field, the fermion number charge must be topological, and 
thus the bosonized fermion number current must be trivially conserved Q We are not 
aware of any theorem that states this, but this is the case in all the known examples 
when either exact bosonization has been performed (1+1 dimensions) [Q 0, or fermionic 
excitations are known to exist in the spectrum of a local bosonic theory in 2+1 [^, or 
3+1 |^T| dimensions. Another observation which indirectly supports this assumption is 
that fermionic fields are not local fields in the usual sense of the word. They satisfy 
local anticommutation, rather than commutation, relations. If these operators are to be 
constructed in terms of local bosonic variables, they must create a nonlocal configuration 
of the bosonic field. If the fermion number bosonizes into a topological charge of the 
bosonic theory, this will necessarily be true. 

Given this assumption, the simplest system that suggests itself as a convenient object 
for bosonization is QED with one two-component Dirac fermion. The reason is that 
Maxwell's equations, 

r = ^d^F"" (4) 

imply that, if one takes as basic variables electric field Ei and the vector potential Ai, 
the fermion number is trivially conserved, due to the antisymmetry of Ff^,y. The fermion 
number charge is topological, since it is equal to the surface integral of Ei at spatial 
infinity. Since there are no global flavor symmetries in the model, it is the simplest of 
its kind. 

In this paper we perform the analog of the Mandelstam construction in the one 
flavor QED3. The paper is organized as follows. In Section 2 we explicitly construct 
the two - component Dirac spinor ipa in the Hamiltonian formalism in terms of bosonic 
variables Ei, and their conjugate momenta Ai. Since the formal definition of ipa suffers 
from the same kind of ambiguities as in the Chern - Simons case, we carefully specify 
the regularization procedure necessary to make it well - defined. In physical terms, the 
fermionic field is an operator which creates an electric charge e at the point x, and also 
a pair consisting of a magnetic vortex and an antivortex of half integer strength vr/e, at 
infinitesimal separation. 

In Section 3 we use the explicit form of the fermionic operators to calculate the 
fermionic bilinears ipTip, and, in particular, show that the operators ipa solve Gauss' 
constraint. The calculation is performed in close analog with the 1+1 dimensional case, 
by expanding all quantities in inverse powers of the ultraviolet cutoff. As in 1+1 di- 
mensions, higher order terms in this expansion contribute finite renormalization of the 

^ In this respect our approach is similar to p9| , although in this reference the resulting bosonic theory 
is nonlocal. 
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coefficients of the lower order terms. As a result, Lorentz invariance has to be invoked in 
order to determine some coefficients in the expressions for the spatial components of the 
current. We then discuss a certain subtlety encountered in our calculation, and which is 
not present in 1+1 dimensions. This problem is due to the presence of the dimensional 
coupling e in 2+1 dimensions, whereas the corresponding couplings in 1+1 dimensions 
are dimensionless. As we show, this restricts the validity of our bosonization procedure 
to scales lower than n, which is related to the UV cutoff A hj fi^ = e^A^. For finite 
this is not a real restriction, since the scale fi becomes infinte in the infinite cutoff limit. 
However, this restriction prevents us from taking the limit at finite A. This we 

believe, however, is not an artifact of our procedure, but rather a necessary consequence 
of the fact that the theory contains only one two - component fermion. The change of 
dynamics at scales higher than /i is necessary to avoid the fermion doubling problem. We 
calculate the algebra of the bilinears, and find that it reproduces the tree level algebra 
with the addition of Schwinger terms. These terms are also seen in perturbation theory 
at the one loop level. 

In Section 4 we find the bosonized expression for the energy - momentum tensor. 
Again, as with the spatial components of the current, we use Lorentz invariance to 
determine several coefficients. It is shown that, with the right choice of these coefficients, 
the resulting theory is Lorentz - invariant in the limit A ^ oo, and one recovers Maxwell's 
equations as its equations of motion. 

Section 5 is devoted to the discussion, and some extensions of our results. In partic- 



ular, we indicate how the well known regularization ambiquity of one flavor QED3 p2 
appears in our framework. It is possible to construct a modified Fermi field Xa, which 
solves the modified Gauss' constraint x^X = ^/^diEi + en/2TTB. The coefficient n must 
be an integer, to ensure the correct quantization of the fermion number charge. The new 
bosonized theory differs from the old one precisely by an addition of a Chern-Simons 



term to the Hamiltonian 23 . 



In Appendices we discuss the simple intuitive mechanism underlying the anticommu- 
tation relations of ip, discuss the geometrical aspects of our construction, and give some 
details of the calculation of the bilinears. 



2 The construction of Fermi operators 

The 2+1 dimensional quantum electrodynamics is defined by the following Hamiltonian 
(we work in the timehke, or Weyl, gauge Aq = 0), 

H = ^E^ + ^B^ + tpYiidi + eAi)ij + m^^ (5) 
together with the Gauss' constraint, 

d,E, = ei,^ (6) 
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The two component Fermi field ipa, ot = 1,2, and the bosonic variables Ei, and Ai satisfy 
the canonical (anti) commutation relations, 

{Mx),My)} = SapS\x - y); [E\x),A\y)] = t6'^\x - y) (7) 

As is usually the case with gauge theories, the Hamiltonian eq. acts on a large Hilbert 
space, which contains unphysical degrees of freedom. Those must be ultimately elimi- 
nated by solving the constraint eq.(^). One usually solves the constraint by expressing 
the longitudinal part of the electric field in terms of the matter fields. The complete set of 
operators that are then used to span the physical Hilbert space consists of the transverse 
part of Ei, and the matter fields ip. In the context of our problem, however, we view 
eq.(^) as one of the bosonization equations. We will therefore retain both components 
of the electric field, and instead solve the Gauss' constraint by constructing the doublet 
of anticommuting fermionic operators ipa in terms of Ei, and Ai |. 

Substituting those back into the Hamiltonian eq.®, we will obtain the completely 
bosonized form of the theory, defined on the physical Hilbert space of QED3. Our 
bosonization procedure is therefore defined only for gauge invariant quantities. We will 
also calculate the fermionic bilinears, and find that they satisfy the tree level algebra, 
modified by Schwinger terms. The appearance of these Schwinger terms is also seen in 
perturbation theory at the one loop level. 

To construct the operators ipa, we must first fix the gauge freedom associated with the 
time - independent gauge transformations, generated by the Gauss' constraint. We do 
this by considering ip in the Coulomb gauge. Those are the gauge - invariant operators. 



where 



ie / d?yei{y - x)Ai{y) 



(8) 



^^(^) = (9) 



^_Xi 

is the electric field of a point charge. The exponential factor ensures the gauge invari- 
ance of the operator ■0*^'^. In the following we will omit the superscript CG, and always 
understand that the fermionic operators we are constructing are defined by eq.(^. It is 
important to note, that the Hamiltonian eq.(|^) can be rewritten in terms of these oper- 
ators. In this formulation, the covariant derivative in eq. contains only the transverse 
part of the vector potential Ai, Aj . 

In addition to solving the constraint eq.(|^), the fermionic operators ipa must satisfy 
the following conditions [|: 

i. Carry unit electric charge, [4'aix),diEi{y)] = etlJa{x)S'^{x — y); 

^ Since the fermion field in QED3 carries only one global quantum number, namely, the electric charge, 
it should be representable in terms of one bosonic field. This bosonic field will not be free of course, 
but will rather be strongly interacting. Our construction gives precisely this counting: the fields Ei and 
their conjugate momenta describe a photon, and an additional scalar degree of freedom. 

■^Bosonization of the massive Schwinger model in 1+1 dimensions can be formulated in precisely the 
same terms. 



6 



ii. Transform correctly under rotations. For convenience, we choose the basis of Dirac 
matrices in which the rotation generator is diagonal, 
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(10) 



Then, 



^i-^e'^/Vi; ^2-^e-''^/V2 (11) 
where is the rotation angle; 

iii. Fermionic bilinears must be local operators, [^/'^■?/'/3(a;), 0(?/)] = for x 7^ and 
for any local gauge invariant operator 0{x). 

To satisfy the first condition, we take the following ansatz. 



Ipaix) = kAVa{x)^{x)Uaix) 



where, 



$(x) 



exp 



ie / d'^yciiy - x)Ai{y) 



(12) 



(13) 



Here A is the ultraviolet cutoff, /c is a finite dimensionless constant, which depends on 
the precise definition of the UV cutoff, and the operators Va and Ua commute with the 
charge density operator. To ensure the anticommutativity of ip^s, we take the following 
forms for f/„ and Kv, 



Vi(x) = —iexp 
Uiix) = 



1- J d'y{0{x - y) ~ ir)dMy) 
-^J d'yeiy - x)miy) 



V,{x) = -tV^x); 
U2{x) = Ul{x) 



(14) 



where 9{x) is the polar angle at the point x. 
With these definitions we find, 



'^i{x)i'i{y) 

il)i{x)i)\{y) 
'ip2{x)i)2{y) 
i^2{x)'^l{y) 

i^i{x)^2{y) 
iJi{x)iiJ2{y) 



7/;i(?/)V^i(x)e^'('(^-")-''("-J^)) = 
V^l(y)V^i(x)e*(^(^"^)"^(^"J'» = - 
^2(y)V^2(x)e'('(^-")-'("-'^)) = - 
4(y)V^2(x)e~'^('(^-")-'("~'^» = 

4(y)^i(a;)e"'" = -ijl{y)Mx] 
^2iy)i>iix)e''' = -ip2iy)^iix) 



-^i(l/)^i(a;) 
iljl{y)^i{x); 
'^2{y)ip2{x) 
-■^2{y)'^2{x) 



(15) 



The factor kA in eq.(|T2D ensures the correct dimensionality of the fermionic fields, and 
the correct normalization of the anticommutators at coincident points. 



{4{x), My)} = {4{x), My)} = s\x- y) 



(16) 



Here a comment is in order about the precise meaning of eqs.(|l^-[T^). The function 
6[x) is defined only modulo 27r. On the other hand, the eigenvalues of the charge density 
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operator are quantized in units of 1/e. As a consequence, the operators Ua and Va (and 
thus ipa) are double - valued. This, in fact, is to be expected on the general grounds. 
Under rotation by 271 the fermionic operators change sign, but the bosonic operators Ai 
and Ei are unchanged. Therefore ip cannot be a single valued function of A^, Ei, and the 
coordinate x P3 . Further geometrical details of the construction, and its interpretation 



along the lines of are contained in the Appendix B. 

These manipulations are, however, too formal, since the commutators of V{x) and 



U{x) with <l>(x) in eq.(|T2|) are singular. The singularity is due to the fact that ^{x) 
creates a charge at the point x, and consequently the commutator involves the ill - 
defined factor e*^(°)/2^ The expression for the fermionic operators needs regularization. 
The natural way to regularize it is by point splitting. We will therefore consider the 
following expressions, 

iP'Kx) = Vi{x + 7])^{x)Ui{x - r]); (17) 

1p2{^) = V2{x + ri)^{x)U2{x — 1]) 

The length of the regulator t] is taken to be proportional to the inverse of the UV cutoff 
|?7| oc 1/A. However, considering rj in any fixed direction breaks the rotational symmetry 
of the problem. To restore rotational invariance, we average over the direction of 77. The 
averaging should be performed with an appropriate phase factor, in order to ensure the 
correct rotational properties of the operators. Under rotation by an angle 0, 6 transforms 
as follows, 6{x) —* 6{x) + 0. We arrive therefore at, 

ipiix) = lim — / dfie~^~^ijji{x), ^'2(2;) = lim — / dfje'^'^^jj^i^) (18) 

A-»oo 2lT J A->oo 2tT J 

where the integral is over the angle of the vector 77, and r) is the unit vector parallel to 
1]. It can be explicitly checked, that, 

{M^),My)} = 0^ \x-y\»l/A (19) 

Therefore, in the limit A 00, we regain the standard anticommutation relations. 

However, we are still not done with the construction of tp. The point that remains 
to be settled is the following. As defined up to now, the fermionic operators depend 
only on the longitudinal component of the electric field E^". Any bilinear that one 
would calculate using these expressions would also depend only on Ej^, and not on the 
transverse part Ef. However, E^" is not a local field, and it is therefore impossible to 
obtain local bilinears with this definition of fermionic operators. (We have checked this 
by explicit calculation.) It turns out, that this problem can be remedied by modifiyng the 
expressions for V{x) and U{x), so that they become creation and annihilation operators 
of a magnetic vortex of half - integer strength, 

Vi{x) = -^exp 1^ 1 d^y [{e{x - y) ~ 7T)d^E,{y) + 27TG^^\y - x)e,,d,E,{y)] | ; (20) 
f/i(x) = exp I d^y [eiy - x)d,E,iy) + 2nG^^\y - x)e,,d,E^iy) 
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V2{x) = -iv,\xy, U2{x) = ulix) 

Here G^'^\x — y) = — ^ \n(jj^x^), x"^ = x^x\ is the Green's function of the two - dimen- 
sional Laplacian, with IR cutoff /i. 

The physical meaning of the operators and Ua is clear from the following commu- 
tation relations ||25|| , 



[V^ix),Biy)] = —Viix)5''ix-yy, [U,ix),Biy)] = -U,{x)5\x - y) (21) 
e e 

This modification does not change either the anticommutation relations or the rotational 
properties of ip. Moreover, now the point splitting procedure used in regularizing if) has 
a very natural interpretation, since the vortex operators V and U are the dual variables 
of QED^ . In any lattice regularized version of the theory the charged field $ should 
live on the lattice sites, and the vortex fields should live on the sites of the dual lattice, 
and therefore be point split from $. 

Equations (ID, (|18D and are our final expressions for the fermionic operators in 
terms of the bose fields Ej and Aj. 



3 Calculation of bilinears 

Our next step is to calculate fermionic bilinears which do not contain derivatives. As 
usual, those should be defined with the help of a point - splitting procedure. We use the 
following definition, 

Mx) = i!(x)r^(x) = — I rfeV^r(^~) 

Svr J 

(22) 

It is implicitly understood in eg. (|2^) that the limit A — oo is taken at the end, after 
(independent) averaging over the directions of e and rj. The insertion of Wilson factor is 
appropriate for the definition of bilinears in a gauge theory. Since we are constructing 
the Coulomb gauge fermions, it is only the transverse part of the vector potential that 
appears in the Wilson factor. The phase expjzxr} is inserted to project onto the relevant 
irreducible representation of the 2D rotation group, while averaging over e. Thus, for 
r = 7o and F = 1 we have x(e) = 0, while for F = 7+ = 71 + Z72, x(f) = ^(^)) and for 
F = 7_ = 7i — ^72, x(f) = •~^(^)- The ratio rj/e is arbitrary, but the final results should 
not depend on it. 

The calculation of bilinears proceeds in complete analogy to the H-l dimensional 
case. Namely, we expand the expression eq.(p2l) in powers of the inverse cutoff, and 
retain only the terms that do not vanish in the continuum limit A 00. This procedure 
has a certain caveat that one has to keep in mind. The operators which are multiplied 
by inverse powers of A, and are formally small, may, in fact, give finite contribution 
in the continuum limit, if these operators have high enough dimensions. This problem, 
in fact, arises also in 1-1-1 dimensions, where the effect of the higher order terms is to 



^/'^(x + e),7°FV'(x-e) 



ie I dxiA^ 
I Jx-e » \ 
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renormalize the coefficients of the lower dimensional operators in a way consistent with 
the symmetries of the problem 0] . We will return to this point later, and argue that the 
situation in our case is similar. 

Let us begin by calculating the charge density and the mass operators. For that we 
need iplipi, and 

where the c - number phase -Ri is given by, 

Ri{e, e, V) = - 2e) - + 2e) + 9{r] - 2e) - 9{r] + 2e) (24) 
and the functions appearing in eq.(PB|) are defined as, 

My) = ei{y + e) - ei{y - e); (25) 
aliy) = e{y-e + O-d{^ + ^-y) + 0{v-e-y)-e{y + e + r]) 
^aUy) ^ G^'\y - 6 + - G^'\e + ^ - y) + G^'\v - e - y) - G^'\y + e + 

We now have to expand the operatorial part in the exponential in eq.(p3D in Taylor 
series in e^, ^j, and rji. The only subtlety here is that the derivatives do not commute 
when acting on 6{x), and thus the order of derivatives has to be specified. Keeping in 
mind the physical picture, however, it is clear that we must ffist expand all expressions 
in powers of rji and C,i at fixed e^, and only afterwards expand in powers of e. The order 
of taking the derivatives is therefore unambigous. Expanding eq. (|23|) up to terms of the 
second order, we obtain (for details see Appendix C), 



(26) 



Jii = ViVi = 27rA;2A2 
Here Ei = €ijEj, and the averaging is defined as follows, 

<a>=^J dedr/deae-^[''('')-^«)llme-i«i(^'«''^) (27) 

Performing the averaging, and choosing the constant k appropriately, we get, 

^iVi = j^diEi ~A-E (28) 

The calculation of ip2'^2 proceeds analogously, and gives the parity conjugate of eq.(p^), 

4'4'2 = ^diEi + A-E (29) 
2e 

We obtain therefore for the charge density Jo and the mass term J, 

Jo = V-V = -^^Ei (30) 
e 
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J = iIjiIj = -2A ■ E (31) 

The correction to these expressions are of the order 1/A^, and we neglect them. 

We now turn to the calculation of the spatial components of the current. Following 
the same steps as in the previous derivation, we find, 



(32) 



The phase R2 is given by, 

^2(e, V) = 0{i - 2e) - e{i + 2e) - ^(r^ - 2e) + e{r, + 2e) (33) 

and, 

ol{y) = 9{y-t + i)-e{e + i-y)-e{'n-e-y) + e{y + t + 7^) (34) 
^al{y) ^ G^^\y - 6 + - G^'\t + i - y) - G^'\v - e - y) + G^'\y + e + 

This again is to be expanded in powers of the inverse cutoff. After some algebra (see 
Appendix C), and keeping two terms in the expansion, we obtain. 



J- 



21 



-/c^A^2 < e,- > eA- - 



(35) 



e 4e^ 

- < eiejek > didjAk ^ < eiejek > AiAjAk 



1 ^ ^ i ^ - ■ 



where the averages are now defined as. 



< a >= 



(36) 



Already at this stage we see that the calculation of the spatial components of the current 
is more involved. The leading term is proportional to the UV cutoff. The next - to - 
leading term in the expansion is of order 1/A, but its commutator with the leading 
term is finite. We therefore have to keep at least those next order terms which give a 
contribution to this commutator, even though they vanish in the naive continuum limit. 

We thus have to come to grips here with the problem mentioned at the beginning of 
this section, namely, that discarding the formally small higher order terms in the above 
expansions is too naive. The simplest manifestation of this problem is the following. 
Explicitly evaluating the averages in eq.(^), and keeping only the leading terms, one 
obtains. 



Ji = ip-fitp = -enAAi H (3EiiAE) + 

eA L 



(37) 



However, the constants k, j3 and 7 in the above formula depend on the ratio of the 
regulators |?7|/|e|. On the other hand, it is clear that the final result may not depend 
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on this ratio. This means, that the formally small terms that we discarded do not 
disappear completely, but renormalize the constants in eq. (|37D . The situation here is 
precisely analogous to the 1+1 dimensional case There, the same naive expansion 
procedure gave an incorrect overall scale of ji, which resulted in jo and ji not forming 
a Lorentz vector. In 1+1 dimensions one could, in principle, explicitly evaluate all the 
corrections coming from the higher order terms. Instead, the requirement of Lorentz 
invariance was sufficient to fix the scale of ji. In the present case we use precisely the 
same method. In the next section, we will calculate the energy - momentum tensor, and 
require that the theory be Lorentz invariant in the continuum limit A ^ oo. We will 
also require that the spatial components of the current, together with the mass term, 
satisfy (up to possible Schwinger terms) the tree - level algebra. 

Now it remains to convince ourselves, that the contribution of the higher order terms 
leads to finite renormalization of the coefficients in the expressions for the Jj's. In 1+1 
dimensions this is a trivial consequence of the fact that the coupling constant in, say 
the sine - Gordon theory, is dimensionless. Since the coupling constant e in QED3 is 
dimensionfull, we cannot use similar arguments here. In fact, at first glance, it seems 
that some of the terms we have discarded are actually more important than the ones we 
have kept. The expansion in powers of 1/A will bring down terms of the same general 
form as in eq. (P7[) , but multiplied by powers of the factors of the following three types. 



A 



eA 



(38) 



The first type of factor is harmless, since, clearly, it can lead only to a finite renormal- 
ization of lower - dimensional terms. Since the scaling dimension of the vector potential 
in perturbation theory is 1/2, the second factor is small at any physical scale, and can 
be discarded. The third term, however, looks dangerous. The perturbative dimension of 
Ei is 3/2, and this factor seems therefore to be of order A^/^. If that were true, it would 
imply that the terms containing higher powers of this factor are more important than 
the lower order termsQ. 

This conclusion is not correct, however, for the following, rather subtle, reason. It 
turns out that, apart from the UV cutoff A, the bosonized theory has an additional UV 
scale /i = (e^A^)^/^. The appearance of this scale can be seen as follows. The fermionic 



operator defined in eq. (|TB|) is essentially a product of a vortex and an antivortex operator, 
separated by A^^. The UV scaling behavior of the vortex operator in perturbation theory 
is known P3|. At short distances it scales as an exponential. 



< V{x)V*{y) >oc \x — yl^exp 



— X 



y\ 



(39) 



where c is a constant, and a/2 is the scaling dimension of the electric field, < E{x)E{y) >(x 
^ Assuming the perturbative scaling of the electric field at short distances, a = 3, 



\x-y\' 



* Although in QED2 the couphng constant has the dimension of mass, the electric field does not 
have a transverse component. Therefore the scaling dimension of E/e is zero, and the problem does not 
arise. 
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we find that the fermion operator that we have constructed, at distances larger than 
1/A, but still small relative to any physical distance scale, behaves, roughly, as. 



< ^lix)^dy) >^ 1^ - y\ 



exp 



e'^A?\x — y\ 



The scale /i is therefore clearly a crossover scale. At distances larger than. 



1 

7' 



(40) 



(41) 



the exponential in eg . (|40|) can be expanded in power series. In that case the fermion 
propagator scales as a power, which is consistent with its perturbative behavior. 

At distances smaller than l//i the non - point - likeness of ip becomes important. 
As a result the scaling behavior at these short distances in the bosonized theory must 
be different from the perturbative one. In fact, if the fermionic operator is still to scale 
with a power law for << 1/x « A, the leading UV behavior of the propagator 
of the electric field must be < Ei{x)Ei{y) >|^._j^|<i/^~ \x-y\^ • ^^^^ return to this 
point later, and see that this short distance behavior emerges naturally from the bosonic 
Hamiltonian. 

Assuming for the moment that this is indeed the correct asymptotics, we find that 
our procedure of calculation of the bilinears is indeed self - consistent. The order of 
magnitude of the "dangerous" correction factor is. 



eA eA 



0(11 



(42) 



and the corrections due to the higher order terms in the expansion of the bilinears can 
only lead to a finite renormalization of the lower order terms. 

We therefore take the current in the general form. 



Ji = ^p-fiij = -enKAi + — : (3Ei{AE) + -fE^A, 

eA L 



(43) 



Clearly, to make this expression well defined, we must normal order it with respect to 
the perturbative vacuum. This implies the subtraction of all terms with at least one 
Wick contracted pair of operators. It has been rigorously shown that, in 1+1 and 2+1 
dimensions. Wick ordering with respect to the free measure of fields in an interacting 
theory always removes the most divergent parts of the correlation functions We 
assume that, in the present case, this is sufficient to subtract the divergent part of the 
operator in the square brackets in eg. (|43|) . 

To be more specific, we assume the following form of the UV asymptotics of the corre- 
lators of the bosonic fields, which is the most general one consistent with the UV scaling 
dimensions discussed above, rotational symmetry, and parity transformation properties. 



lim^_y < Ai{x)Aj{y) > 



Ti- 



\x - y\ 



+ 2r; 



{x-y)i{x - y)j 
\x — y\-^ 



(44) 



13 



|x — \x — 



\x — y\^ \x — y\ 



en „ (X — y)i(x — y)j — (x — y)i(x — y)i 

\x — y\ \x — y[^ 

The expectation values of the quadratic operators are then given by, 

^ < ^2 >= 2(gi + g2)A'; < >= 2{r^ + r^)^ (45) 

< A- E>= 2(si + S2)A^; < A- E >= -2m{pi + p2)A 
< >= -(ri + 2r2)A^ 

The numerical coefficients are therefore subject to the following conditions, 

gi + 92 > 0; ri + r2 > 0; n + 2r2 < 0; si + = (46) 

The first three conditions are obvious, while the last one should hold in a Lorentz in- 
variant theory. It follows from the fact that E ■ J is the zeroth component of a Lorentz 
vector, and therefore must have a vanishing expectation value. 

In order to fix the coefficients in eq.(^), we ffist impose the condition that, together 
with other normal ordered bilinears, this current satisfies the tree - level current algebra, 
up to possible Schwinger terms, and terms that vanish in the continuum limit, 

[Ji{x), J,{y)] = '2te,,J{x)6\x - y); [J,(x), J{y)] = -2te,,Jj{x)6\x - y) (47) 

Evidently, with the mass bilinear given by eq.(0), the second equation holds without 
any corrections. Calculating the ffist commutator, we obtain, 

[J,(x), J,(y)] =2e,,{-2(gi+g2)(pi+P2)(/3+27)'+[-/t(3/3+27) + (gi+g2)(/3+27)'] : A-E : 

(48) 

+ ^^[-m(pi + P2)(/? + 27) : E' : +^ -.E^A-E-]} 
This gives two conditions, 

/3 + 27 = 0; (5=- (49) 
The current is then determined as, 

J, = ^jiij = -enKAi + ^ \2Ei{AE) - E^A^ + 4(pi + p2)mK^ (50) 
The constant n will be determined in the following section. 
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The fermionic bilinears are now expressed as local functions of Ei and Ai. Their 
algebra can be calculated explicitly, 



[Ji{x),Jj{y)] =2ieijJ{x)5'^{x -y); [Ji{x) , J (y)] = -2ieijJj{x)5'^{x - y) (51) 



[Jo{x), Ji{y)] = -i 



2i 



d]6\x 



[J(x), My)] = -Ei{x)e,,d]6\x - y) 
e 



The first two commutators are the canonical tree - level ones. The other two exhibit 
explicitly the Schwinger terms mentioned above. This, again, is a common feature of 
our procedure and the bosonization in 1+1 dimensions: the Schwinger terms, which in 
perturbation theory appear at the one loop level, appear in the bosonized theory in the 
tree level commutators. 



4 The Hamiltonian and Lorentz invariance 

We now calculate the energy - momentum tensor of the theory in terms of bosonic fields. 
The gauge invariant, symmetric energy-momentum tensor of QED3 is given by, 

with the bosonic and fermionic parts, 

T^^ = ^{iPYD''^ + i^l^'Di'ilj - D^ijjYi^ - D^'tPYi^) (54) 
Here D = d — leA is the covariant derivative. In order to bosonize the fermionic part 



of the energy - momentum tensor, one has to calculate the bilinears in eg. (|5^) . Their 
direct calculation in terms of the fermionic operators constructed in Section 2 gives the 
general form of the different terms which enter the Hamiltonian, but cannot establish 
their coefficients (because of the unknown finite renormalization due to the higher order 
terms). 

In 1+1 dimensions there is an alternative way of calculating T^*^, using its represen- 
tation in the Sugawara form. The energy - momentum tensor in a free fermionic theory 



can be written as a suitably regularized product of the currents p8|. In 1+1 dimensions 
it is, 

r ~\ IT) — 

T^^ = [{ JM, r} _ ^MVV,] + g^^-i^l^, ^} (55) 

Thus the knowledge of the bosonized form of the currents suffices to obtain the bosonized 
energy - momentum tensor in this case. 
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In 2+1 dimensions one can still write a formal analog of eq.(E 



"A 



1 r 



{r, r} - g^^j'A- : { r} ■. +9''' ■. J'Jx ■] + g'-'ji^P, ^} (56) 



Here : S : means normal ordering with respect to the perturbative vacuum. The normal 
ordered term is formally of order 1/A, and vanishes in the continuum limit. However, 
being a composite operator of high dimension, it cannot be discarded a priori. The 
simplification in 1+1 dimensions is that the normal ordered part vanishes, due to the 



symmetries of the 1+1 dimensional Dirac matrices PS|. This does not happen in higher 
dimensions, and we cannot use the Sugawara construction directly. 

We can still use the construction, however, as a guide to the form of T^'^. The 
naive (i.e. neglecting the normal - ordered contributions) Sugawara form of T'^'^ contains 
precisely the same terms that one gets from the explicit calculation. We will therefore 
take an ansatz which is consistent with both approaches up to finite renormalizations, 
and determine the coefficients by requiring that Lorentz invariance is recovered in the 
continuum limit. 

The general form for the Hamiltonian density is, 
T°o = + Ie' + ^ {d,E,)' + b^AA' + 1 [c : (A ■ E)2 : +rf : E^A' ■]+fA-E (57) 

If the theory is to be Lorentz invariant with standard transformation properties for 
Ei,B,jQ, and Jj, Maxwell's equations should be satisfied in the continuum limit. We 
therefore require, 

B{x) = t[H, B{x)] = -ei.d^Ej + 0(1/A) (58) 

Ei = i[H, E,] = -Ji + tijd.B + 0(1/A2) (59) 

Since Jj itself contains terms of order 1/A, we demand that the second equation be 
satisfied to order 0(1/A^). 

Calculating the commutators, we find the following conditions on the coefficients 

/ = 0, 6 = ^, c=--, d=^ (60) 
2 K 2k 

This still leaves the coefficients /t and a undetermined. Our next step will be therefore 
to impose the complete Poincare algebra. 

A sufficient condition for Lorentz invariance of a theory invariant under spatial rota- 



tions and translations is the following commutation relation |29 



-z[Too(x),Too(2/)] = -{T^{x)+T^{yW5\x-y) (61) 



^In fact, for these values of the coefficients the inhomogeneous Maxwell's equation is obtained exactly, 
without corrections 0(1/A^). The correction term to the homogeneous equation is of the form, 0(1/ A) — 
+2 / nke^^ dk[{A^ E^) + {AE)A- — {MKk/2)A-)]. Normal - ordering the operator in the square brackets, 
we find that its divergent part vanishes. Thus, even for finite A, this equation retains the form of a 
conservation law. 
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It can be easily verified (by multiplying eg. (|6TD by linear functions of x and y and inte- 



grating) that eq. (|6lD implies the closure of the Poincare algebra [29|, with the following 
genertors: the energy / (i^xT°°(x), the momentum Pj = / (PxTf, and the angular mo- 
mentum and boost generators, given, respectively by. 



Li — j d x^XqT^ XjT L — j d xa^jX^Tj 



(62) 



It should also be verified that the angular momentum defined in eq.(|62D does indeed 
generate rotations when acting on Ei and Ai. Calculating the commutator in eq. (|6lD, 
we find for the momentum density. 



T° = BEi - 2aKdiEiA. 



(63) 



4a(gi + q2) + ri + 2r2 ~ 4a(gi + gs) - ri - 2r2 
GijA ■ E OijA ■ E 



K 



+^^L±^a,(A,i, + A,A,) + 0(1/A) 

In arriving at this expression we have extracted finite parts of the formally vanishing op- 
erators using the operator product expansion (normal ordering) with the UV asymptotics 
of the propagators given in eq. (^5]) . The relevant Wick contractions are, 



< diEi{x)Ej{y) > Oiy) 



lim 

lim^_y^o < diEi{x)Aj{y) > 0{y) 
lim^_y^Q < ekidkAi{x)Aj{y) > 0{y) 
lim^-y_,o < ekidkAi{x)Ej{y) > 0{y) 



{q. + q^Wdpiy) 
1 



(64) 



(Pi + P2)mAejidiO{y) 
Vi + 2r2)AejidiO{y) 



^{Pi + P2)mAdjO{y) 



Requiring that the angular momentum generator has the standard "orbital" and 
"spin" parts. 



L 



d^xeijXiEkdjAk — A ■ E 



we obtain from eq. (|63|) the following relations. 



1 



+ ^2) 



or, usmg eq. 



2fi;' ri + 2r2 

A < ^2 > 



< 52 > 



(65) 



(66) 



(67) 



The expression for the momentum density is also consistent with what is expected in 
a theory of a local vector field A^. The momentum operator Pi = J d^xT^'^[x) generates 
spatial translations by the usual rule. 



[Aix),P,\ 



-idiAj{x) 



(6^ 
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(and the same holds for Ei). 

We have now determined all the coefficients in the Hamiltonian density and the 
current in terms of ^§2^ and < A ■ E >. Our final expressions are (we record their non 
- normal - ordered forms) , 

= '-B^ + '-E^ + ^{d.E^r + "^^AA^ (69) 



^ -2{A- Ef + E^A^ +MA-E; 



2kA 

T^i = BE, - id,Ej)A, + ^diiAE) + ^djiAAj + AjA.,) + O(^); (70) 

o A 



Jj = —enAAi + 



en A 
with. 



2EiA ■ E - AiE' 



M ~ 
+ —Ei 
e 



M = -^<A.E>=^^^^m (71) 
kA k ^ ' 

A-^K 

and K defined in eg. (|67D . 

We want to stress here that, although eq. (|69D contains a parameter n, it does not 



define a one - parameter set of theories. The bosonized version of QED3 corresponds 
to a unique choice of n, which satisfies eq.(|67D. Unfortunately, since eq.(|69D defines a 
strongly interacting theory, we cannot determine the numerical value of n. This would 
involve the solution of the model (at least in the UV region) for arbitrary n, calculating 
< and < E'^ >fj, and solving the selfconsistency equation (|57|). 

We end this section with a comment. In arriving at the form of eq. (|57D , and subse- 
quently eq.(^), we have truncated the series for the Hamiltonian density at the order 
1/A. This amounts to neglecting terms of the form ^2x3-^^^^- This is consistent with the 
asymptotic form of the correlation functions given in eq.(|45|), since the power counting 
based on it tells us that we have kept all the relevant terms. We can estimate the scale 
at which the irrelevant terms become important. Assuming perturbative scaling at low 
energies, a simple scaling argument shows that these terms are unimportant relative to 
the terms we kept in eq.(p7D at energy scales smaller than fj. = (e^A^)^/^. This is again 
consistent with our previous estimate (see Section 3). At this scale all terms become 
equally important, and, once this happens, the scaling dimensions of fields will cease to 
be canonical. It is generally the case that the magnitudes of irrelevant and relevant terms 
in a Hamiltonian coincide at the UV cutoff scale. In our case this translates into the 
statement that, near the scale A, the electric field must have dimension 1. Although this 
argument does not constitute a proof of eq.(^5D (and the proof cannot be given without 
solving the theory), it shows the selfconsistency of our assumptions. 
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5 Discussion 



In this paper we have performed the analog of Mandelstam's construction in 2+1 dimen- 
sional QED3, that is, we have constructed the two - component Dirac spinor field entirely 
in terms of the bosonic fields Ei, and their conjugate momenta Ai. The fermionic bilin- 
ears are local functions of the bosonic variables. They satisfy a current algebra which 
includes Schwinger terms. The bosonized theory is Lorentz invariant in the continuum 
limit. These aspects of our construction are very similar to what one encounters in 1+1 
dimensions. There are, however, certain features which are conceptually different, and 
which reflect the differences between 1+1 - and 2+1 - dimensional physics. We now make 
several comments on these issues. 

Firstly, the fermionic operators eq.(0) anticommute only at distances larger than the 
ultraviolet cutoff. In fact, the phase in the anticommutation relations contains a power 
tail of the form (l/(A|a; — y\)^, where s is a number of order one. This is closely related 
to the fact (discussed in Sections 3 and 4), that the construction involves two UV scales 
A and fi oc (e^A^)^/^, and that the equivalence with continuum QED^ holds only below 
the lower scale yU. This also means that, if we were to discretize our bosonic theory, 
the fermionic operators would not have canonical local anticommutation relations. The 
lattice spacing a should be identified with the smallest distance scale required in our 
regularization procedure, that is, a oc 1/A. The Fermi fields would then anticommute 
only at distances larger than = a. This is contrary to the assumptions made in the 
proof of the Nielsen - Ninomiya theorem and provides a way in which the theory 
avoids the standard fermion - doubling problem. In fact, recall, that a Hamiltonian 
theory of one staggered lattice fermion field with local anticommutation relations in 2+1 
dimensions becomes in the continuum limit a theory of four Fermi degrees of freedom. 
Our theory, however, has only two fermionic degrees of freedom in the continuum limit. 
This aspect is different from the 1+1 dimensional case, where a theory of one staggered 
lattice fermion leads in the continuum limit to a theory of one Dirac fermion. Accordingly, 
the anticommutation relations in 1+1 dimensions are canonical at all distances up to the 
cutoff. 

Secondly, we comment on the explicit appearance of the UV cutoff in our final ex- 
pressions, which is also a novel feature. The UV cutoff never appeared explicitly in 1+1 
dimensional formulae, since there the theory of free fermions bosonized onto a theory 
of a free bosons. In higher dimensions there is however no reason to expect that the 
same will happen. On the contrary, one expects that fermions free in the UV will be 
represented by strongly interacting bosons. Indeed, the bosonic Hamiltonian eq.(|69D 
contains, apart from quadratic terms, also a quartic interaction term with a coefficient 
of order 1/A. Since this coefficient scales with the inverse power of the cutoff, by naive 
power counting, the term is irrelevant. However, this term is very important in deter- 
mining the scaling dimensions of various operators. For example, if we were to omit it, 
the Hamiltonian would describe a (non - Lorentz - invariant) theory of two free bosonic 
degrees of freedom A^, with scaling dimensions one. Both components of the electric field 
would then have dimension 3/2. We know, however, that the divergence of electric field 
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(which is proportional to the fermion number density) must have scahng dimension 2. 
The formally irrelevant quartic interaction term is precisely responsible for the required 
change of the scaling. 

It is also clear that this term cannot be treated in perturbation theory, since, even 
though it appears to be irrelevant, its effect is not small. The situation here is very 
similar to the one in some strongly interacting 2+1 - dimensional theories defined at 
finite fixed points of the coupling constant. The examples that recently have attracted 
much attention are the four - Fermi theories The four - Fermi interaction term is 



perturbatively irrelevant in 2+1 dimensions. However, if the coupling constant is defined 
at (or infinitesimally close to) its fixed point value, a renormalizable and nontrivial theory 
is obtained. The scaling dimensions of various operators at the fixed point are different 
from their values in the free theory. Since the coupling constant has the dimension of 
inverse mass, its fixed point value is always of the form g = j^, where A is the UV cutoff, 
and is a pure number. The theory eq.(|5UD should be understood in the same sense. 
It is defined at the fixed point value of the coupling of the naively nonrenormalizable 
quartic interaction term, which is responsible for the correct scaling of the fields. In this 
spirit eq.(^), although it was not derived here from corresponding (3 - function, should 
be understood as a fixed point condition. The appearance of the positive powers of A 
in eq.(|69|) is also natural, once one realizes that the theory is interacting. Generically, 
in an interacting theory any coupling constant with positive dimension must scale as 
a power of the UV cutoff. This is the origin of the cutoff in the term in eq. (|69|) . 
The appearance of the cutoff in the expression for the conserved current eg . (|3^) is then 
inevitable, since the form of the currents is determined by the Hamiltonian. 

Unfortunately, the very fact that the bosonic theory we obtained is strongly interact- 
ing prevents us from analysing these questions quantitatively. Approximation schemes 
that were useful in analysing the four - Fermi interactions, such as the 1/N expansion, 
or the e - expansion, cannot be applied here, since an extension of the model either to 
large number of fermionic species, or to dimensions different from 2+1 is not available 
at the moment. 

Even so, there are several directions in which our approach can be extended. The first, 
and most straightforward one, is to understand the regular izat ion ambiguity of QED^ 



2^ . The fermionic operators we have constructed solve the Gauss' constraint. If one 
relaxes this condition, there are additional possibilities p3|- natural modification 
is to substitute for Ei in eq.(pUD the linear combination Hi = Ei + ne'^eijAj. This gives 
X^X = \diEi + KeB. (x is the modified Fermi field). However, the coefficient n is not 
totally arbitrary in this case. It is crucial for the derivation that / (Pxdjli has quantized 
eigenvalues. Otherwise, fermionic operators will not be double valued, and the bilinears 
will not be local bosonic operators (see Appendix C). Since the magnetic flux in QED^ 
is quantized in units of the coefficient k must be given by with n - an integer. 
This possible modification is the reflection of the well known regularization ambiguity 
in fermionic QED3 which leads to the appearance of the induced Chern - Simons 
term in the action or, in the Hamiltonian formalism, to the modification of the Gauss' 
law constraint. This modification of fermionic operators will also modify the expressions 
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for the fermionic bilinears, and the energy - momentum tensor. The resulting bosonic 
theory can be obtained by the methods used in this paper. 

Thirdly, although the theory we have derived is Lorentz invariant, it is not written 
explicitly in terms of Lorentz covariant fields. It would be desirable to make one more 
step, and find a formulation in terms of a different, covariant bosonic field. We believe 
that the most natural candidate for this is the magnetic vortex field. In 2+1 dimensions 
it is a scalar field It has the additional virtue that in terms of it not only is 

the electric current trivially conserved, but also the electric charge has a meaning of 
a topological winding number, and is quantized classically The vortex operators 

appeared naturally in our construction of the Fermi fields. The counting of the number 
of degrees of freedom suggests that this additional transformation should be possible, 
since the complex vortex field is equivalent to two real bosonic fields. It would also 
be interesting to extend the construction to theories with more fermionic species. This 
generalization may be simpler in terms of the scalar vortex variables. 

Finally, we note that the fermionic operators which we have constructed involve a 
product of vortex and antivortex operators, carrying opposite magnetic fluxes. Therefore 
the operators do not carry net magnetic flux, and thus are not analogous to 2 + 1 
- dimensional dyons. The mechanism of their anticommutativity is different from that 
discussed in , and subsequently extensively exploited in the analysis of Chern - Simons 
theories. In fact, the approach to the problem of bosonization presented here has a 
natural generalization to 3+1 dimensions. The fermionic operators can be constructed 
quite easily in a form analogous to eq.(|18]). The anticommutation relations can be 
achieved by the mechanism described in the Appendix A. While in 2+1 dimensions the 
fermionic operator creates a point charge and an infinitesimally close pair consisting of 
a magnetic vortex and an antivortex of half - integer strength, in 3+1 dimensions it 
creates a point charge and an infinitesimally small vortex loop, again of half - integer 
strength. By averaging over the orientations of the loop, and by ensuring the correct 
transformation properties under the axial transformations, it is possible to construct the 
four component Dirac spinor, and calculate the bilinears. The work along these lines is 
currently in progress . 
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6 Appendix A. A simple mechanism for anticommu- 
tation. 

In this Appendix we give an intuitive explanation of the basic element in the construction 
of the anticommuting operators ipa, and point out its similarity to the 1+1 dimensional 
case. Recall that, in 1+1 dimensional Abelian gauge theory, the fermionic operators can 
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be expressed in terms of the gauge field as follows, 



4'i,2{x) = exp 



—le 



dyA{y) 



exp 



'^i-E{x) 
e 



(73) 



The two factors in eg. ([73|) have a very simple meaning. The first exponential factor 
creates a unit charge at the point x, 

^9E(z)e''/-oc'^^^(2') = \dE{z) + 5(x - z) (74) 



e e 
The second exponential (up to a surface term) can be rewritten as 



V± = exp 



-oo T 
J2 J- 



(75) 



It therefore measures the total electric charge in the half space to the right of the point 
X. To see why two such operators anticommute, consider ipa{x)ipi3{y). Since either x > y 
or y > X, only one of these two operators creates the electric charge in the region where 
the other one is measuring it. When we change the order of the fermionic operators, 
the expression picks a phase according to e^e^ = e^e^e^^'^\ But always, only one of 
the V's measures the charge. As a result the phase is always ±7r, and the operators 
anticommute. 

It is clear therefore that in any number of dimensions one can achieve anticommuta- 
tion relations (but not rotational invariance) from a construction of the following kind. 



'ip±{x) = exp 



i / d'^yei{x-y)Ai{y) 



exp 



d^y-miy) 

M e 



(76) 



where ei{x) is the field of a point charge {diei{x) = 6'^{x)), and the volume M over which 
the integral is performed in the second factor is the half space "to the right of" the point 
X, defined for example by yi > xi. The objects constructed in this way anticommute 
precisely for the same reason as in 1+1 dimensions. 

Although 2+1 dimensional fermionic operators which we discussed in this paper where 
not represented in the form of eg. (|TBD , the reason for their anticommutativity is essen- 
tially the same. Forgetting momentarily about the regularization, and the factor ordering 
(which are crucial for the calculation of bilinears, but not for the anticommutation prop- 
erties of the operators at distant points), we can write eg.(p!7D as. 



ipa{x) oc exp 



i / d'^ye,{x - y)Ai{y) 







exp 





d^y{d{x - y) - d{y - x) 



^))miy) 
(77) 

For our present simplified discussion, we can think about the polar angle 6{x) as a 
function with a cut. The cut should begin at the origin, and go to spatial infinity. Let 
us choose the second axis as the direction of this line. Then, for [x — y)i > 0, we have 
9{x — y) — 0{y — x) = —TT, and for (x — y)i < we have, 9{x — y) — 0{y — x) = n. The 
second phase factor in eg.([77|) then becomes precisely the integral over the half space 



(x — y)i > of the electric charge density \diEi, with the coefficient vr. 

We therefore see, that the anticommutativity in our construction is achieved by pre- 
cisely the same mechanism as in 1+1 dimensions. 
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7 Appendix B. The geometry of our construction. 



In this appendix we discuss the geometry of our construction of the fermionic operators, 
and point out its connection to the formulation of Finkelstein and Rubinstein |^ . 

The function 9{x) is defined by the Cauchy - Riemann equation [0, 



e,,d,G^^\x) = ^d,e{x) (78) 

Here G^'^\x — y) = — ^ ln(/i^a;^). The differential equation is solved by, 

9{x) = I dyie,^% (79) 
Jc{M,x) y' 

where the curve C{M,x) starts at a base point M, and ends at the point x. It is 
convenient to take the base point M to infinity. The function 6 then depends only on 
the point x, and the first homotopy class of C, i. e., the number of times the curve 
C winds around 0. A point on the universal cover can be parametrized by x, and the 
homotopy class of this curve. Therefore ^ is a single valued function on the universal 
covering space Moreover, e*2, and therefore also ipaix), is a single valued function 



on the double cover of the plane R^. 

We wish, however, to consider our theory as QED^ on a plane, rather than on its 
cover. There is a natural mathematical framework - introduced by Finkelstein and Rubin- 
stein - which, suitably generalized, accomodates this interpretation of our construction. 

To set the scene, we will now discuss the relevant setup in some generality. 

Let X be the (time 0) configuraton space of the theory. We will assume, that X 
is a Banach space, noncompact in the natural topology. Let the set {0i}i=i,... comprise 
of quantum fields, i.e., the operator - valued distributions on X. We shall assume that 
there is an action of the full Lorentz group on the fields (pi, which are bosonic quantities. 
The Fock space over X, , is the field manifold of the theory. Let Q denote the set 
of all quantum fields on X. Q is explicitly realized as the Li^ space of the functional 
Gaussian measure, formally, 

AT n n (^H^) [- < 0; K<p >] (80) 

zex i 

for the positive definite kernel K defined by the quadratic part of the action. The 
integration is over S {X), which is the space of Schwartz distributions on X. 

Consider now a multivalued functional of the fields 0, say ^ will in general 
depend not only on the fields themselves, but also on continuous functions {Oj}j=i,...N^ 
defined on the universal cover X of the configuration space. We shall, however, elect 
to work not with the cover, but with X itself, accordingly, we will take the domain 
of \1/ to be Q X n~^{(C(X)^}, where is the natural projection from the cover to 
the configuration space. In this approach, the "functions" 9 are multivalued. We shall 
denote the domain of the functionals \l/ by Q. 
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Now, let Fa be a flow on Q. Any such flow induces a continuous 1 - parameter 
family of linear transformations on the functionals. We will take here, for definiteness, 
the flow to be induced by rotations about a fixed axis. Accordingly, we shall take = 1, 
and 9i = 6, where ^ is a polar variable parametrizing a 1 - dimensional submanifold, 
perpendicular to the rotation axis. The linear transformation is then exhibited as, 

^{4>,e)^^{F-'4>F^,e + a) (81) 

We shall take a e [0, 27r], accordingly Fq — 1, is a rotation by 27r, and ^2^0 = 
0. However, is assumed to be multivalued, which implies, that 7ri(Q) ^ 0. Since 
Q is a rather complicated infinite - dimensional set, we usually cannot compute this 
homotopy group directly. However, as we will show, we can relate this group to a much 
simpler homotopy group. To this end, consider a given topological sector of the theory, 
characterized by the asymptotic behavior of the operators, 

M^) ^ <Poi (82) 

Denote the corresponding Fock space superselection sector by J-'^{<^q). Now, consider a 
loop in Q. This can be viewed as a mapping of the unit cell J^^icpo), with 

values of the fields and functions on 91'^^"^^+^ specified. Hence, 

{^""M) (83) 

as is easily seen by considering = 0o- 

It is now obvious, how this construction should proceed in our formulation of QED^. 
We take X — H'^, X — C/Z2, and = {E,A). Since we are working in the gauge - fixed 
formulation, the corresponding Gaussian measures are well - defined. We have explicitly 
construced a pair of double - valued functionals ^, which change their signs under a 27r 
rotation. We thus deduce, that 7Ti{Q) 7^ 0. This implies, by the isomorphism shown 
above, that 7r3(jF"^(0o)) does not vanish. We conjecture, that this result is tied to the 
fact, that while the operators Ei have continuous spectra, the operator diEi has a purely 
discrete spectrum. This imphes, that a "naive" picture of the resulting Fock space is 
certainly not correct, and this space has nontrivial topological characteristics. 



8 Appendix C. Calculation of fermionic bilinears. 



In this appendix we give details of the calculation of fermionic bilinears. We start with 



l{x + e),ip — l{x — e) 



I^M?l,kl«i/A 



(84) 



We have to represent the product of fermionic operators as a single exponential, expand 
it in powers of the regulators, and average over the directions of the regulators. Using the 



24 



Baker - Campbell - Hausdorff formula, e^e^ = e^+^+ 2 ["^'^1 , which is valid when i^] 
is a c - number, we find, 

= ^--mn)-m-\R^i^,i,v\^F,iA,E) (gg) 

Here the c - number phase Ri is given by, 

Ri{e, C, V) = die - 2e) - 9iC + 2e) + ^(77 - 2e) - 9(7) + 2e) (87) 
and the operatorial part is, 

Fr{A, E) = eJ d^yh{y)A,{y) + d^yal{y)d,E, - a}r{y)t,,d,E,{y) (88) 

with, 

f,{y) = ei{y + e)-ei{y-e)- (89) 
ai(y) = %-e + 0-^(e + e-y) + ^(r?-e-y)-% + e + r7) 

i-a^(l/) ^ G^^\y - e + - G(^)(e + ^ - 1/) + G^^^r^ - e - y) - G^^\y + 6 + 77) 

We now expand the operatorial part in Taylor series. Since the derivatives do not com- 
mute when acting on 9{x)^ one has to specify their order. However, physically it is clear 
that we must first expand in 77 and ^ at fixed e, and therefore the order of derivatives is, 
in fact, unambiguously determined. We use the following identities, 

J d'yejiy - x)Aj{y) = -Af{x) (90) 
d: J d'y9{y - x)d,E,{y) = 27r^^ 
df ld'yG'{y-x)ejkdjE,{y) = 

where the superscripts L and T denote the longitudinal and transverse parts of the fields, 
respectively. Expanding the operatorial part up to the terms of second order, we obtain, 

^iF,iA,E) ^ ^ ^ 2ieeiAf - i—{r] - OiEi - i—iv + Oi^A^i - 2ehiejAfApl) 

- 47re,{7j - OjA^Ej - ^(r; - Oiiv - OAEj 
^iej:jy.Ajiy) ^ I + 2ieeiAj - 2ehiejAj AJ 

The phase Ri depends on the ratio |77|/|e|. It turns out, however, that changing this 

iRl 

ratio results in multiplying Ime 4 by a constant. Since this is the only quantity that 
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appears in the present calculation, its variation can always be compensated by a suitable 
redefinition of k in eq. ([T8D . We therefore evaluate Ri for |e| >> \ri\, 



Ri{e,v,O = 



eijei7]j > 0, eijei^j > 

eijeir]j > 0, eijei^j < 

tijeifij < 0, eijti^j < 

eijeir]j < 0, eijei^j > 



(92) 



The next step is to calculate the averages, defined as. 



< a >= 



in-" 



(93) 



We find, 



< 1 >= 0; < r), ± >= 0; < (r), ± ii){f]j ± Q >= 0; < e, >= 0; < e.ej >= (94) 



Assembling the terms, we obtain. 



(95) 



-i<{ri + ^),ej > -djEi- < ei{rij - ^j) > {Ai, Ej} 

= ^A'\e\\v\[ld.E,-2A.E] 

The calculation of ip24'2 proceeds in the same fashion. The only difference is that the 
electric field enters all expressions with the opposite sign, and the averaging is done with 
the complex conjugate phase factor. 



< a> -z 



in-" 



As a result the averages now become. 



We therefore get. 



1 - i 



i;l^2 = ^A''\eM[-diEi + 2A-E] 
2 e 



Choosing the length of the regulators in an appropriate way, 

\e\\r]\ = k-^A-^ 

we obtain, 

ip^^ = -diE„ t/jtp = -2A ■ E 
e 



(96) 

(97) 

(98) 

(99) 
(100) 
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We now turn to the calculation of spatial components of the current, 



J_(x) = 'ipl{x)'ijji{x) 



1 



in 



i)l{x + e), - e) 



ie r^~^^ dxiA^ 

5 Jx—e * 



Following analogous steps, we find. 

The operatorial phase is, 

F2{A, E) = iej £yf,{y)A{y) -^J <fyal{y)d,E, + al{y)e,,diEj{y) 

and, 

al{y) = e{y-t + i)-e{t + i-y)-e{'q-t-y) + e{y + t + 7]) 



l'?Mei,kl<xi/A 
(101) 



(102) 
(103) 

(104) 



(105) 



—al {y) ^ G(^) (2/ - e + - G^^^ {e + i-y)- G^'^ {v - e - y) + G^'^ (y + e + v) 
Zn 



The c - number phase R2 is given by, 

R2{e, e, V) = - 2e) - + 2e) - ^(r^ - 2e) + ^^(r/ + 2e) 



^2(e,r/,O = 
i?2(e,r7,0 = -27r 

^2(e,r/,O = 
/22(e,r7,0=27r 



eijei?7j > 0, eijCi^j > 

UjUVj > 0, eij-ej^j < 

eijeiT]j < 0, eij-ej^j < 

eijeir]j < 0, eijei^j > 



(106) 



(107) 



The Taylor series of F2 starts with a term of order zero, 

- - / d'y[e{y) - e{-y)]d,E,{y) 
e J 

Again, using the Baker - Campbell - Hausdorff formula, we obtain, 

giFa ^ _giF2 + i/ d''yl9{y)-ei-ymEM^-lJ d?y[e(y)-e{-y)]d,E,(y) 



:io8) 



(109) 



The function 6{ii) — 0{—y) is equal either to tt or to — tt. The last operatorial factor in 
eq.( |109| ) can therefore be written as. 



exp 



d^y[e{y) - d{-y)]d,E,{y) 



exp[-m{Q+ - Q-)] 



(110) 



where (5+ is the total electric charge in the regions of space where 9{y) — 9{—y) = tt, 
and Q_ the electric charge in the regions where 6{y) — 6{—y) = —tt. Multiplied by the 
factor exp{— ivrQ} in eqs. (|102|-|103D, the operatorial factor becomes. 



exp{— i27r(5+} 



:iii) 
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Since the electric charge is quantized in sets which are sufficiently regular, so that their 
Newtonian capacity can be defined, this operator is, in fact, the identity. Note that the 
quantization of charge is crucial for the locality of J,, since otherwise the exponential in 
eq. (|lll|) would be a complicated nonlocal operator. Now we can perform a straightfor- 
ward Taylor expansion of the operatorial phase in eqs.( |102[ - |103D , and average over the 
directions of the regulators. The nonvanishing contributions are. 



j_ = <ti> eA. 

4e3 



(112) 



1 ^ ^ i ^ ^ ■ 

- < eii^j + rij){^k + rik) > AiEjEk + ^ < (^ii^j - Vj){^k + Vk) > diEjEk 



where the averages are now defined as. 



< a >= 



(113) 



Calculation of J+ proceeds in an analogous fashion. The expression for the current 
obtained in this way is, 

Ji = i)-^ii) = eAk'^A^\t\Ai+ (114) 



B{djdjAi + 2didjAj) + Ce^AiA^\ + —K^\tM^ \DEi{AE) + TE'^Ai 



The coefficients A, B, C, V and can be calculated explicitly. We do not give their 
values, however, since they are regularization - dependent. The final expression for Jj 
is determined in Sections 3 and 4 from the requirement of Lorentz invariance and tree - 
level current algebra. 
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